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We investigate the zero-temperature behavior of the classical Heisenberg model on the triangular 
lattice in which the competition between exchange interactions of different orders favors a relative 
angle between neighboring spins $ G (0, 27r/3). In this situation, the ground states are noncoplanar 
and have an infinite discrete degeneracy. In the generic case, i.e. when $ 7^ 7r/2, arccos(— 1/3), the 
ground state manifold is in one to one correspondence (up to a global rotation) with the set of 
non-crossing loop coverings of the three equivalent honeycomb sublattices into which the bonds of 
the triangular lattice can be partitioned. This allows one to identify the order parameter space as an 
infinite Cayley tree with coordination number 3. Building on the duality between a similar loop model 
and the ferromagnetic 0(3) model on the honeycomb lattice, we argue that a typical ground state 
should have long-range order in terms of spin orientation. This conclusion is further supported by 
the comparison with the four-state antiferromagnetic Potts model [describing the $ = arccos(— 1/3) 
case] , which at zero temperature is critical and in terms of the solid-on-solid representation is located 
exactly at the point of roughening transition. At $ 7^ arccos(— 1/3) an additional constraint appears, 
whose presence drives the system into an ordered phase (unless $ = 7r/2, when another constraint 
is removed and the model becomes trivially exactly solvable). 

PACS numbers: 75.10.Hk, 75.50.Ee, 05.50.4-q 



I. INTRODUCTION 

The recent discovery of the compound NiGaS4, and 
the suggestion that it might be a spin liquid or a spin 
nematic,^ has revived the interest in the Heisenberg 
model with bilinear and biquadratic interactions defined 
by the Hamiltonian 

H - ^ JiS, • s,, + J2 (s, • s,o' , (1) 

to') 

where {jj') stands for nearest neighbor pairs on a lattice. 
In one dimension, the spin-1 case has been very thor- 
oughly investigated following Haldane's prediction that 
the pure Heisenberg model is gapped,^ and a rather com- 
plete picture of its properties has emerged both in zero 
(see e.g. Refs. [3-5] and references therein) and in finite 
magnetic field. 

In two dimensions, most of the attention has also been 
devoted to the spin-1 case, and some trends have emerged 
from the investigation of this model on the triangu- 
lar and square lattices. With the standard notation 
J I — J cos 9, J2 = J sin 9, the phase diagram as a func- 
tion of 9 consists of four main phases: an antiferromag- 
netic phase (with two sublattices on the square lattice 
and three sublattices on the triangular lattice) around 
9 = 0, followed counter-clockwise by a three-sublattice 
antiferroquadrupolar phase up to = 7r/2, a ferromag- 
netic phase around 9 = t:, and finally a ferroquadrupolar 
phase that persists until the antiferromagnetic order sets 
in. 

In comparison, little attention has been paid to the 
classical limit, where spins are classical vectors of length 



1. The main reason lies in the fact that, unlike for spin- 
1/2 models, in the spin-1 case the classical limit is not 
equivalent to the Hartree approximation (in which the 
ground state is approximated by a product of local wave 
functions) since local wave functions are not necessarily 
purely magnetic but can also describe quadrupolar states. 

In the classical limit, the most unusual feature in- 
troduced by the presence of the biquadratic exchange 
is the possibility (realized in the parameter range 
— 2J2 < Ji < 2 J2 with J2 > 0) to have the minimum of 
the interaction of two spins in a noncollinear configura- 
tion in which they make an angle 

$ = arccos (^~^^ ^ (0' (2) 

with respect to each other, while the overall orientation 
of the spins is arbitrary.^^'^^ Such an interaction may be 
called a rotationally invariant .spin- canting interaction. It 
should not be confused with the Dzyaloshinskii-Moriya 
interaction, which also induces a canting between a pair 
of spins but forces them to lie in a specific plane. 

Note that the physics will be essentially the same for 
any rotationally invariant canting interaction between 
pairs of spins that leads to the same angle $, for in- 
stance interactions that include higher powers of Sj ■ Sj', 
in which case the relation (2) would have to be replaced 
by a more complex one. For this reason we prefer below 
to characterize the interaction in terms of the angle $ 
and not in terms of the parameter 9. The classical XY 
model with rotationally invariant spin-canting interac- 
tion on the square lattice has been investigated by Lee et 
al.^^ and the classical Heisenberg model on the triangular 
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lattice by Kawamura and Yamamoto.^^ 

According to Ref. 15, the zero-temperature phase di- 
agram of the classical version of model 1 on the tri- 
angular lattice includes four phases (see Fig. 1): a fer- 
romagnetic phase for tt — arctan(l/2) < 9 < 37r/2, 
a nematic-like phase with coUinear orientation of spins 
for —tt/2 < < — arctan(2/9), a three-sublattice an- 
tiferromagnetic phase for — arctan(2/9) < 9 < 7r/4, 
and a phase with non-coplanar orientations of spins 
for 7r/4 < < TT — arctan(l/2). This phase diagram is 
reminiscent of the spin-1 phase diagram, but the non- 
degenerate ferro- and antiferroquadrupolar phases of the 
spin-1 case are replaced by highly degenerate magnetic 
phases. The nematic-like phase is thoroughly discussed 
in Ref. 15, with the conclusion that it has the same de- 
generacy and correlations as the antiferromagnetic Ising 
model on the triangular lattice. 

The discussion of the non-coplanar phase is much more 
sketchy however. The authors point out that there is a 
degeneracy due to the two possible chiralities of a lo- 
cal umbrella [formed on each triangular plaquette when 
$ G (0, 27r/3)], but no attempt has been made to iden- 
tify the structure of the ground state manifold or to study 
the consequences of the degeneracy on the spin correla- 
tions in a typical ground state. In the present paper, we 
construct the complete classification of the ground states 
which allows us to conclude that the non-coplanar phase 
is a non-trivial one even at zero temperature because the 
total number of ground states grows exponentially with 
the area of the system, as in various versions of the an- 
tiferromagnetic (AF) Potts model and other frustrated 
models with a finite residual entropy per site.^^"^^ In sit- 
uations like this, it is not evident a priori whether the 
degeneracy leads at zero temperature to the disordering 
of the system or to an algebraic decay of correlations, or 
whether long-range order is nevertheless present. 

Interestingly enough, in the two-dimensional antiferro- 
magnets known to have an extensive residual entropy, all 




e = - 31/2 



FIG. 1: Phase diagram of the classical bilinear-biquadratic 
Heisenberg model on the triangular lattice. 



three scenarios are realized. In particular, both the AF 
Ising model on the kagome lattice^® and the three-state 
AF Potts model on the square lattice^^ at zero tempera- 
ture are disordered, the AF Ising model on the triangu- 
lar lattice^^ and the three-state AF Potts model on the 
kagome lattice^*^'^" are critical, whereas the three-state 
AF Potts model on the dice lattice has a genuine long- 
range order. The main aim of our paper is to establish 
which of the three scenarios is realized in the present 
case. 

The paper is organized as follows. Sec. II is devoted 
to a complete classification of the ground states of the 
model in terms of the nonintersecting closed loops living 
on the bonds of the triangular lattice, which allows us to 
redefine the problem as a loop model. In Sec. Ill we show 
that the order-parameter manifold of this loop model has 
the same topology as that of the 0{n) model with n ~ 3, 
which suggests that they belong to the same universal- 
ity class and, therefore, that both of them are in the or- 
dered phase. In Sec. IV, the same conclusion is confirmed 
by establishing a relation with the four-state AF Potts 
model on the same lattice and by relying on some known 
properties of this Potts model. For completeness. Sec. 
V briefly discusses the case $ = 7r/2, where the model 
is trivially exactly solvable, and Sec. VI summarizes the 
results. Appendix A is devoted to the analysis of a sim- 
plified problem corresponding to freezing the spins on 
one of the three sublattices in a perfectly ordered state. 
It also includes estimates from below and from above of 
the residual entropy of the full problem. 



II. GROUND STATE CLASSIFICATION AND 
LOOP MODEL 

We assume that classical spins Sj (three-dimensional 
vectors of unit length) are located at the sites j of the 
triangular lattice and that the energy of the interaction 
of two neighboring spins E{S,S') depends only on the 
angle between S and S' and is minimal when this angle 
is equal to $ G (0, 27r/3). In terms of the coupling con- 
stants Ji and J2 introduced in Eq. (1) this corresponds 
to — 2J2 < Ji < J2 (see Eq. 2). For such values of it 
turns out to be possible to minimize E(Sj,Sji) simul- 
taneously for all pairs of neighboring spins {jj'), which 
means that the system is not frustrated. However, as we 
shall see, it is very highly degenerate. 

The simplest ground states of the model can be con- 
structed by choosing three unit vectors S'^, and S*^ 
which form angles $ with respect to each other. After 
that, one can partition the triangular lattice into three 
triangular sublattices (denoted below A, B, and C) and 
set all spins on sublattice A to S^, on sublattice B to 
and on the third one to S'-^. In such a way, one ob- 
tains a periodic ground state which has a three-sublattice 
structure. Below we call such ground states regular 
states. 

The family of regular states is characterized by an 
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SO {3) X Z2 degeneracy, where the group 5*0(3) is related 
to the possibiUty of simultaneously rotating all spins and 
the group Z2 to the possibility of choosing the sign of the 
mixed product of the three spins located in the corners of 
a given plaquette to be either positive or negative. How- 
ever, it is clear that the family of ground states is much 
wider than the set of regular states, and that, in addition 
to these two symmetry-induced degeneracies, it possesses 
also an infinite set of discrete degeneracies not related to 
global symmetries. Indeed, for any pair of spins S"* and 
forming an angle $ € (0, 27r/3) with respect to each 
other there are always two possibilities to choose the ori- 
entation of the third spin in such a way that it forms 
the same angle $ both with S"^ and with S^. These two 
orientations are mirror images of each other with respect 
to the plane formed by S"* and S^. Since in a regular 
state all plaquettes contain the same triad of spins (S"^, 
and S*^), it is clear that in such a state one can al- 
ways flip an arbitrary spin into its mirror image without 
increasing the energy of the system. 

Flipping one spin leads to a new ground state in which 
six neighboring triangular plaquettes (surrounding the 
flipped spin) form a domain of a different regular state. In 
order to obtain a larger domain, one has to flip other spins 
belonging to the same sublattice. As soon as a domain of 
another regular state is sufficiently large, a new domain 
can be formed inside it by fiipping spins belonging to 
another sublattice, and so on. 

A convenient method of describing the full set of 
ground states is based on the notion of a zero-energy 
domain wall.^^ In the system under consideration, a zero- 
energy domain wall can be defined as a line passing 
through the bonds of the triangular lattice in such a 
way that each segment of this line separates two plaque- 
ttes corresponding to different regular states. This means 
that each segment of such a wall separates two unparallel 
spins which on the original triangular lattice are next-to- 
nearest neighbors and therefore belong to the same sub- 
lattice. Then it is clear that the spin configuration in a 
ground state is uniquely defined as soon as one specifies 
the orientations of three spins on an arbitrary plaquette 
and the positions of all zero-energy domain walls. For 
brevity, in the following we simply call these objects do- 
main walls. 

By definition, a regular state contains no domain walls, 
whereas flipping of a single spin leads to the formation 
of the simplest closed domain wall consisting of six seg- 
ments joining each other at angle 2tt/3, i.e. a hexagon. 
By flipping a larger number of spins belonging to the 
same sublattice one can construct more complex domain 
walls. However, it is clear that any wall has to be closed 
(or end at the boundary) and its neighboring segments 
always have to form angles of ±27r/3 with respect to each 
other, as in Fig. 2(b). 

Since each domain wall can be associated with spin- 
flipping on a particular sublattice of the triangular lat- 
tice, all domain walls can be partitioned into three classes 
(which below are called colors). The domain walls of a 



given color live on one of the three equivalent honey- 
comb sublattices into which the bonds of the triangular 
lattice can be partitioned. Therefore, just by construction 
walls of the same color have no possibility of crossing or 
touching each other. 

A less evident property is that for a generic value of 
$ walls of different colors are also unable to cross each 
other. This is so because when going around an inter- 
section of two walls [see Fig. 2(e)] one has to perform 
successive spin flips on sublattice a (where a = A, B , C) , 
then on another sublattice /3 and once again on a and 
then on /3, and return after that to the same configura- 
tion. At a formal level, this corresponds to the fulfillment 
of the condition 



(3) 



where Ra denotes the operation of spin-flipping on sub- 
lattice a {R^ = 1). However, since each of the operators 
Ra and Rp involves a rotation of a spin around the third 
spin of the triad by the angle 25'($), where 



arccos 



cos $ 
1 + cos $ 



e (tt/B,^) 



(4) 



the application of {RpRa)^ rotates the whole configura- 
tion by 4^'($). It is easy to check that one returns to the 
initial configuration only for $ = tt/2 (when is also 
equal to 7r/2), whereas for all other values of $ the oper- 
ator {RjjRa)^ brings a triad of spins into a configuration 
rotated with respect to the original one, as illustrated in 

Condition (3) can be also rewritten as RaRp = RpRa, 
which means that the zero-energy domain walls can cross 
each other only when operators R^ commute with each 
other. However, walls of different colors can touch each 
other, as shown in Fig. 2(c). 

As a consequence of the rules discussed above, for each 
site of the triangular lattice there are only three possibil- 
ities: (i) the absence of any domain wall passing through 



(a) 



(b) 



(c) • 

yf >j< 




(e) 



FIG. 2: (Color online) The loop representation: (a)-(c) al- 
lowed configurations of loops, (d) additional vertex allowed 
only at $ = arccos(— 1/3), (e) the vertex corresponding to an 
intersection of loops (allowed only at $ = 7r/2). The domain 
walls on different sublattices are shown with lines of different 
types. 
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FIG. 3: (Color online) From left to right, configurations ob- 
tained after successive flips of spins A, B, A, B. In the final 
configuration, the triad of spins is rotated around C by an 
angle 4^($) with respect to the original one. 



this site [see Fig. 2(a)], (ii) the presence of a single do- 
main wall whose segments form an angle 27r/3 with each 
other [see Fig. 2(b)] and (iii) the presence of two domain 
walls touching each other [see Fig. 2(c)]. Therefore, to un- 
derstand the zero-temperature properties of the system, 
one has to study a loop model on the triangular lattice in 
which on each site only the configurations shown in Figs. 
2(a)-(c) are allowed (as well of course as those that can 
be obtained from them by rotation), whereas all other 
configurations are prohibited [for example, the intersec- 
tion of two domain walls shown in Fig. 2(e)]. The two 
exceptions where other configurations of zero-energy do- 
main walls are allowed are the cases $ = arccos(— 1/3) 
and $ = 7r/2, discussed respectively in Sec. IV and Sec. 
V. 

One can introduce on each triangular plaquette k a bi- 
nary variable Xfe = il defined in such a way that it is 
perfectly ordered in any regular ground state and changes 
its sign on crossing each domain wall. In terms of the orig- 
inal spin variables Sj the value of Xk is given by the sign 
of the mixed product Sj^(fe) • [Sj^(fc) x Sj^(k)], where the 
three sites jA{k), jsik) and jc{k) belonging to plaquette 
k are labelled after the sublattice to which they belong. 
Below for simplicity we call variable Xk chirality of pla- 
quette k, although more accurately it should be called 
staggered chirality. 



top panel, the domains corresponding to the same reg- 
ular state are marked by the same number (the num- 
bering being arbitrary). On the bottom panel, the same 
set of numbers is used to mark the nodes of the Cayley 
tree. The distance on the tree corresponds to the minimal 
number of spin flips one has to make (in other terms, the 
minimal number of domain walls one has to cross) to get 
from one regular state to another. 

The tree-like topology of the order-parameter mani- 
fold (that is, the absence of closed loops) is ensured by 
the impossibility (for a generic value of $) to return to 
the same state after crossing a sequence of domain walls 
corresponding to a product of operators Ra which cannot 
be reduced to unity by application of the identity — 1. 
This property follows from the fact that intersections of 
domain walls arc not allowed. 

Another model with the same structure of the order- 
parameter manifold is obtained when one considers a loop 
model on the honeycomb lattice^'^'^'* in which each closed 
loop of length L is ascribed a factor 



w{L) = nK^ 



(5) 



In this model the loops cannot intersect or touch each 




III. ORDER-PARAMETER MANIFOLD 

Let us define neighboring regular states as regular 
states which can be transformed into each other by flip- 
ping the spins of one of the three sublattices. From the 
definition of a domain wall it is clear that in any ground 
state the domains lying on both sides of any domain wall 
belong to neighboring regular states. Since each regular 
state has exactly three neighboring states, at zero tem- 
perature the order-parameter manifold of our system has 
the topology of an infinite Cayley tree (also known as 
Bethe lattice) with coordination number n = 3. Each 
node of this tree represents a regular ground state, while 
the links connect neighboring states. If one ascribes to 
each link on the tree the color of the corresponding do- 
main wall, the three links connecting any node with its 
neighbors will all be of different colors. 

Fig. 4 illustrates the correspondence between the do- 
mains of various regular states and the nodes of the Cay- 
ley tree whose links are colored in such a way. On the 




FIG. 4: (Color online) Top panel: Example of loop covering. 
The numbers refer to the positions of different domains on 
the Cayley tree shown in the lower panel. The numbering is 
arbitrary and has been included to keep track of the relative 
positions on the tree. 
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other simply because the geometry of the honeycomb lat- 
tice does not allow for this. The number n can then be 
interpreted as the number of different colors the loops can 
have. This allows one to interpret these colored loops as 
domain walls^'^ separating the states on a Cayley tree 
with coordination number n by following the same rule 
that the states separated by a domain wall of a given 
color correspond to neighboring nodes connected by a 
link of the same color. 

For n = 3 the only difference between the two models 
is that in our model the loops of different colors are liv- 
ing on three different (interpenetrating) honeycomb lat- 
tices, whereas in the loop model of Ref. 23 they are living 
on the same honeycomb lattice. However in both mod- 
els the loops cannot cross each other or have an over- 
lap, which ensures that the two models have the same 
order-parameter manifold and therefore can be expected 
to belong to the same universality class. 

In Ref. 23 it has been shown that the loop model de- 
fined by Eq. (5) is a dual representation of a ferromag- 
netic 0{n) model whose partition function can be written 
as 



Zo(n) — Tr 



(kk') 



(6) 



where Sfc are n-dimensional unit vectors defined at the 
sites k of the honeycomb lattice, the product is taken 
over all bonds {kk') of this lattice and the trace implies 
the unweighted integration over all variables Sj.. On the 
other hand, it is also well known that in 0{n) models with 
n > 2 the long range order in is always destroyed by 
small continuous fluctuations of the order parameter (the 
spin waves), which leads to an exponential decay of the 
correlation functions (sfc^s/c^) with the distance between 
ki and fc2. 

It follows from the analysis of Refs. 26 and 27 that in 
terms of the loop representation the exponential decay 
of (sfcjSfcj) means that the insertion into the system of 
a string (an open loop) going from fci to ^2 leads to the 
suppression of the partition function of the system by a 
factor which exponentially depends on the distance be- 
tween ki and k2 . This suggests that long loops have large 
effective free energy and accordingly the probability of 
finding a sufRciently large loop also decays exponentially 
with the size of this loop. As a consequence, the typi- 
cal number of irreducible loops separating the two points 
cannot experience an unrestricted growth with the in- 
crease of the distance between these points and has to 
saturate at some finite value. 

This means that in terms of the order parameter de- 
fined on a Cayley tree the system is in the ordered phase, 
that is, it remains localized in some region of the or- 
der parameter manifold even when its size tends to in- 
finity. From the universality we expect the loop model 
constructed in the previous section (which has the same 
symmetry of the order-parameter manifold) to be also in 
the ordered phase. In terms of original spin variables Sj, 



such a situation corresponds to the existence of a genuine 
long-range order. This is especially evident for <I> <C 1. It 
is clear that if in a typical ground state the system is 
localized on the tree within the region of diameter -Dtyp, 
then for $ ^ 1/Z?typ the distant spins will be almost 
parallel to each other. 

To support the hypothesis that there should be no 
qualitative difference between the system in which the 
loops of different colors live on different sublattices and 
the one where they live on the same sublattice, we an- 
alyze in Appendix A a reduced problem with n = 2. In 
particular, we explicitly demonstrate that the system in 
which the non-intersecting loops of two different colors 
live on two interpenetrating honeycomb lattices can be 
in an exact way transformed into the system in which the 
loops of two different colors live on the same honeycomb 
lattice. 



IV. COMPARISON WITH THE FOUR-STATE 
ANTIFERROMAGNETIC POTTS MODEL 

When $ = arccos(— 1/3), different spins Sj can have 
only four different orientations in any ground state, which 
means that at zero temperature the system is equivalent 
to the four-state AF Potts model. It has been shown in 
Refs. 19 and 28 that the set of ground states of the four- 
state AF Potts model on the triangular lattice allows a 
mapping onto the set of ground states of the three-state 
AF Potts model on a fca(?o?7ie lattice, whose exact solution 
at zero temperature was constructed by Baxter in 1970.^^ 
His results imply that the four-state AF Potts model has 
a finite extensive entropy per site whose numerical value 
is approximately equal to 0.379. 

Much later Huse and Rutenberg^" have demonstrated 
that Baxter's solution corresponds to an algebraic decay 
of spin correlations on the same triangular sublattice of 
a kagome lattice: (Sjj^ • 8^2) oc l/r^2 with 77 — 4/3, ri2 
being the distance between ji and j'2 (the application 
of the same approach shows that in terms of the four- 
state AF Potts model on the triangular lattice rj = 1/3). 
This conclusion was reached by constructing another 
mapping^*''^^ which puts each ground state of the three- 
state AF Potts model on the kagome lattice into corre- 
spondence with a ground state of a solid-on-sohd (SOS) 
model describing fluctuations of a two-dimensional in- 
terface in a four-dimensional space. In this SOS model 
height variables Uj (defined on the sites j of the original 
triangular lattice T) are points on an auxiliary triangular 
lattice % defined in the transverse (height) space. The 
only restriction on possible configurations of u is that if j 
and j' are nearest neighbors on T, then Uj and Uji have 
to be nearest neighbors on % . 

To make the description of the system more trans- 
parent, it is convenient to introduce the locally coarse- 



grained heights hfc 



u 



-3A{k) 



U 



Js(fc) 



+ U 



Jc(fe) 



, that is 



the averages of variables u over the three sites belonging 
to a given triangular plaquette of f,'^^'^^ These variables 
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can be considered as defined at the sites k of the hon- 
eycomb lattice H dual to T, and their values belong to 
the honeycomb lattice dual to %. Variables are 
more convenient than variables Uj because, in any regular 
state, all variables are equal (in contrast to variables 
Uj). On the other hand, each domain wall separates 
and hfe' which are nearest neighbors on "H*. Therefore, 
if k and fc' are nearest neighbors on "H, then and h^' 
have to be either equal or nearest neighbors on H*. 

Thus, for $ = arccos(— 1/3), different regular states 
can be associated with different points of 'H^,, which 
means that the order-parameter manifold instead of be- 
ing a Cayley tree with coordination number n = 3 (as for 
generic values of has the topology of a periodic lattice 
with the same coordination number, namely, of the hon- 
eycomb lattice. ^^''^'^ On a formal level, this follows from 
the existence of the identity (RpRa)^ = 1. This difference 
in the topology of the order-parameter manifold leads to 
a change of the universality class, and, instead of being 
long-ranged, the correlations are algebraic. In terms of 
the loop representation introduced in Sec. II, a special 
feature of the case $ = arccos(— 1/3) is that in addition 
to vertices a, b and c (see Fig. 2), the vertex d shown in 
the same figure is also allowed (with the same unitary 
weight as all other vertices). 

It is known from the analysis of Ref. 20 that the 
(2-|-2)-dimensional SOS model corresponding to the 
Potts model under consideration is in the rough phase, 
where the large scale fluctuations of h can be described 
by a free field effective Hamiltonian, 

ffcff = ^ 1 d2r(Vh)2 , (7) 

however, the value of the dimensionless rigidity modulus 
J is such that the perturbation responsible for the dis- 
creteness of h is marginal. In other terms, the model is 
located exactly at the point of the roughening transition. 
Therefore, any additional suppression of height fluctua- 
tions should drive the system into the ordered phase. In 
particular, we argue below that this happens when one 
suppresses the formation of vertices d (which for a generic 
value of $ are not allowed). 

In terms of height variables the main property of the 
loop model defined by vertices a, b and c is that its order- 
parameter manifold does not allow for the existence of 
closed loops in the height space. As a consequence, when 
going around any of these vertices, the vector h makes 
in the height space just a trivial loop that can be con- 
tracted to a point and has zero area. In contrast to that, 
when going around vertex d the vector h sweeps a loop 
whose area is equal to that of the elementary hexagon 
cell of H*. Accordingly, in terms of height variables hfc 
an interaction that keeps the weights of vertices a, b and 
c unchanged but leads to the partial or complete sup- 
pression of vertices d can be written as 

^^ = yEr,^ (8) 
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where > 0, 

1 ^ 

is the area of a loop swept by h when going counterclock- 
wise around the hexagonal plaquette of the dual lattice 
surrounding site j of the original lattice, and p numbers 
(in the same direction) the six sites of the dual lattice 
belonging to this plaquette. Naturally, complete suppres- 
sion corresponds to Jd — -l-oo . 

In the framework of a continuous description, Eq. (8) 
should be replaced by 

where ri and r2 are the two components of the vector 
r and the quantity in brackets is the Jacobian of the 
mapping r — >■ h. It seems natural to expect that the 
addition of such an interaction with > to the large- 
scale Hamiltonian (7) would lead to the suppression of 
the fluctuations of h and indeed the perturbative treat- 
ment demonstrates that already in the first order of the 
expansion in powers of one obtains a correction to J. 
For Jd > this correction is positive and therefore has 
to shift the system into the ordered phase. This provides 
one more argument in favor of the conclusion that for a 
generic value of $ the system has to be in the ordered 
phase, which corresponds to the existence of a genuine 
long-range order in terms of the original spin variables 
S,. 

However, the existence of such a long-range order does 
not necessarily imply the presence of the long range-order 
in chirality. Even if in a typical configuration the system 
is localized in some region of the order-parameter man- 
ifold, its distribution does not have to be centered on a 
particular node of the order-parameter tree, but can also 
be, for example, symmetric with respect to the center of 
a link connecting two neighboring nodes. In such a case 
there will be no long-range order in chirality, although 
the original spin variables Sj will be ordered on all three 
sublattices. 



V. THE CASE $ = tt/2 

In the case <& — 7r/2, the spins on neighboring sites al- 
ways have to be perpendicular to each other (where "al- 
ways" means "in any ground state"). As a consequence, 
each spin always remains flippable independently of the 
flips made by other spins. In such a situation it is natu- 
ral to describe the ground states of the system in terms 
of Ising pseudospins aj — ±1 such that Sj = S^^^^aj, 
where a{j) — A, B,C denotes the sublattice to which 
site j belongs and S'^, S^, S*" are three unit vectors 
perpendicular to each other. Since all possible sets of 
aj are allowed, the residual entropy per site is equal to 
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ln2 « 0.693 and the zero-temperature spin-spin correla- 
tion function (Si • Sj) vanishes as soon as i ^ j. 

However, the system has long-range order in terms of 
the spin-nematic order parameter, S'^S'' — 6°'^/3,^^ with 
three different easy axes (which are perpendicular to each 
other) for the three sublattices. In Refs. 11 and 12 de- 
voted to the investigation of model (1) with 5 = 1, 
the phase with such a structure of the order parameter 
(realized for 7r/4 < 9 < tt/2) was called the antiferro- 
quadrupolar phase. The difference between the quantum 
system with S = 1 and the classical limit, 5 — >■ c», is 
that for S — 1 the order parameter on each sublattice is 
like in an easy-plane spin nematic, whereas for S ^ oo 
it is like in an easy-axis spin nematic. 

As it was already mentioned in Sec. II, in terms of the 
loop representation the case $ = 7r/2 corresponds to the 
situation where the intersections of loops [see Fig. 2(e)] 
are not prohibited. In such a case, the three subsystems of 
loops (on different sublattices) are completely decoupled 
from each other. 



VI. CONCLUSION 

We have investigated the classical Heisenberg model 
on the triangular lattice with a spin-canting interac- 
tion, i. e. an interaction such that the energy of a pair 
of spins is minimized when the angle between them 
is equal to $ e (0,7r). In particular, the interaction of 
two spins has this property if the antiferromagnetic bi- 
quadratic exchange is of comparable strength with the 
bilinear exchange. Our analysis is focused on the case 
$ S (0, 27r/3), where the energy can be minimized simul- 
taneously for all bonds. The family of ground states is 
then characterized by a well-developed degeneracy corre- 
sponding to an extensive residual entropy, which makes 
the question whether at zero temperature the system is 
ordered, disordered or critical a non-trivial one. 

After constructing a complete classification of the 
ground states of the model, we have demonstrated that 
at zero temperature its order-parameter manifold has 
the structure of an infinite Cayley tree with coordina- 
tion number n = 3 and therefore is isomorphic to that of 
the 0{n) model with n = 3. This suggests that the two 
models belong to the same universality class and there- 
fore the system under consideration should be, like the 
0{n = 3) model, in the ordered phase, which means long- 
range order in terms of the original spin variables Sj . This 
conclusion has been confirmed by showing that our model 
can be constructed by imposing an additional restriction 
onto the four-state AF Potts model on the triangular 
lattice, which leads to the suppression of the fluctuations 
of the free field describing the large-scale fluctuations of 
this Potts model and therefore pushes it from the point 
of roughening transition into the ordered phase. 

Numerical simulations of the classical bilinear- 
biquadrartic Heisenberg model on the triangular lattice 
demonstrate^'"' that at J2/J1 = —3 and very low temper- 



atures the spins in a typical configuration form a non- 
coUinear structure analogous to what can be expected 
from a typical ground state in the presence of weak con- 
tinuous fluctuations. However, the authors of Ref. 15 have 
not studied how the spin-spin correlation functions be- 
have in the T — )• limit (at any finite temperature the 
long-range order in spin orientation in any case has to 
be destroyed by the spin waves^^'"^^), so their results can- 
not be used for checking our conclusions. Moreover, it 
still remains to be elucidated if the behavior of the sys- 
tem at very low temperatures is determined mostly by 
its zero-temperature properties, or whether the removal 
of the accidental degeneracy by continuous fluctuations 
(spin waves) can play an important role, as in the case 
of a Heisenberg or XY antiferromagnet on the triangular 
lattice in the presence of an external magnetic field. '^'^ 
Note that the only possibility for the non-coUinear phase 
considered in this paper to have a genuine long-range or- 
der at finite temperatures consists in having a long-range 
order in chirality, and our conclusions allow for the real- 
ization of such a scenario. 

The two exceptions from the generic behavior 
described above are the cases $ = 7r/2 and 
$ = arccos(— 1/3) for which the model at zero temper- 
ature is exactly solvable. For $ = 7r/2 the fluctuations 
of any pair of spins are uncorrelated, (Si • Sj) = 0, 
but the system has long range order in terms of the 
spin-nematic order parameter with three different easy 
axes for the three sublattices (the antiferroquadrupolar 
phase^^'^^). On the other hand, for $ — arccos(— 1/3) the 
model is equivalent to the four-state AF Potts model, in 
which the chirality and spin correlations are algebraic, 
(XfeiXfcfca) <^ l/''i2 (see Ref. 20), {S-j, ■ S^,) cx l/r^f . 
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Appendix A: Reduced problem 

Let us consider a "reduced" loop problem which differs 
from the one formulated in Sec. II by the suppression of 
the formation of loops on one of the three honeycomb 
sublattices. This is equivalent to saying that the spins 
on one of the three triangular sublattices are not allowed 
to flip and remain aligned in the same direction (below 
this sublattice is denoted 7c and two other sublattices 
Ta and Tb). In such a situation, the formation of loops 
on the honeycomb lattice JiAB formed by the sites of 
Ta and Tb is impossible and the loop system consists of 
only two types of loops (on honeycomb lattices Hac and 
-Hbc)- 
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When all spins on 7c point in the same direction, 
= S*^, the orientation of spins on the two other sub- 
lattices (7a and 7b) can be described by a single vari- 
able ipj the angle between the projection of Sj on the 
plane perpendicular to S*^ and some reference direction 
in this plane. Since the angle between neighboring spins 
is always equal to <&, the difference between the values 
of ■0 on neighboring sites of Hab has to be equal (mod- 
ulo 2tt) to ±\E'($), where \E'($) is given by Eq. (4). For 
$ G (0, 27r/3), *($) belongs to the interval (7r/3,7r) and 
is a monotonically increasing function of 

For ^'(<I)) ^ tt/2 (that is, $ ^ 7r/2), the existence of 
such a relation between the orientations of neighboring 
spins allows one to introduce integer variables Uj defined 
on the sites of honeycomb lattice Hab in such a way 
that on neighboring sites of Hab (belonging to differ- 
ent sublattices) they always differ by ±1. In particular, 
one can choose variables Uj to be even on Ta and odd 
on Tb ■ Then on neighboring sites of the same triangular 
sublattice these variables either differ by ±2 (when these 
sites are separated by a domain wall) or are equal to each 
other (in the absence of such a wall). 

Accordingly, the system of AC and BC loops can be 
discussed in terms of a (2 -|- l)-dimensional solid-on-solid 
(SOS) model defined on Hab - In the SOS representation, 
variables Uj have the meaning of surface heights, whereas 
each domain wall plays the role of a step between the 
regions in which the values of either or of dif- 
fer by ±2. The SOS model defined in such a way can 
be considered as a generalization of the BCSOS model 
(introduced by van Beijeren'^^ on the square lattice) to 
the honeycomb lattice. In particular, the representation 
of the model in terms of Ising-type pseudospins aj = ±1 
can be constructed with the help of the same relations as 
for the BCSOS model,^^ 

af = cos[(7r/2)u/] , (Ala) 
af = sin[(7r/2)«f ] , (Alb) 

where we still follow the convention that variables Uj are 
even on Ta and odd on Tb- The pseudospin variables 
CTj are defined in such a way that each flip of the original 
spin Sj leads to a change of sign of aj . Accordingly, pseu- 
dospins of opposite signs belonging to the same sublattice 
are separated by domain walls. 

All configurations allowed in this SOS model enter 
the partition function of the reduced problem with the 
same weight (like in the infinite temperature limit of the 
BCSOS model). The existence of the mapping to a SOS 
model confirms that the suppression of spin flips on one 
of the three sublattices reduces the order parameter man- 
ifold of the system from a tree with coordination number 
n = 3 to a trivial tree with coordination number n — 2. 

The SOS model on Hab introduced above can be fur- 
ther mapped onto the 20-vertex model on the triangular 
lattice introduced by Baxter.'^'' In order to perform the 
mapping, one should put on each bond of Tc an arrow 
directed in such a way that the larger of the two heights 



at the ends of the AB bond crossed by this arrow is al- 
ways to the right of the arrow. Then on each site of Tc 
one will have three incoming and three outcoming arrows 
(the ice rule) which precisely corresponds to the selection 
of the allowed vertices in the 20-vertex model. 

The exact solution of the 20-vertex model on the trian- 
gular lattice found by Baxter"'*'^ demonstrates that when 
all vertices have the same weight (as in our case) the di- 
mensionless free energy per vertex is equal to — ^ In ^ . 
In terms of the current discussion this means that the 
residual entropy of the reduced problem (per site of the 
original triangular lattice T) is equal to | In ^ « 0.318. 
The residual entropy of the full loop problem should be 
larger than this value but lower than that of the four- 
state AF Potts model (approximately equal to 0.379 ^^). 

Let us now compare the contributions to the partition 
function of the reduced problem which can be associated 
with different configurations of the BC loops (in other 
terms, with different allowed configurations of variables 
M^). The contribution from the configuration in which 
tire BC steps are completely absent (that is, all variables 
are equal to each other, = u"^) can be found very 

easily. In such a situation, all variables can acquire 
two values ± 1 independently of each other. 

Therefore, the corresponding contribution is equal to 2^° , 
where Nq is the number of sites on each of the three 
sublattices. 

If there is a closed BC loop of length L separating 
the regions in which the values of the variables differ 
by 2, then on all L/2 sites from Tb passed by this loop 
the variables cannot fluctuate and have to be equal 
to the average between the values of on both sides 
of the loop. This means that the presence of a loop of 
length L decreases the number of allowed conflgurations 
on sublattice Tb from 2^0 to 2^0-^/2-1-1^ rpj^^ additional 
factor of 2 appears because for each loop configuration 
this loop can represent two types of steps (positive and 
negative). 

After noting that an analogous reduction is induced by 
every loop on Hbc, one can conclude that after summa- 
tion over all allowed configurations of variables one 
obtains an SOS model defined on the triangular lattice 
Ta- In this model the values of height variables on neigh- 
boring sites of Ta can either be equal or differ by ±2, 
whereas each step is ascribed a weightfactor exponen- 
tially decaying with its length. This model can also be 
interpreted^'^''^® as the Ashkin- Teller model on the trian- 
gular lattice'^^ with one of the weights equal to zero. 

Thus we have reduced a model in which the loops of 
two different types are living on different interpenetrating 
honeycomb lattices to a model in which all loops are liv- 
ing on the same honeycomb lattice. However the weight 
factor corresponding to each loop 

w{L) ^2{l/^/2)^ . (A2) 

contains a factor 2 related to the existence of two types 
of loops. The loop model defined by the weight fac- 
tor (A2) is just a particular case of the more general 
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model^'^ with 'w{L) given by Eq. (5). For — 2 < n < 2 
this model is known^'* to have a phase transition at 
Kc{n) 1/^2 + 72^. 

This suggests that for n = 2 and K = 1/^2 the SOS 
model on Ta defined two paragraphs above is located 
exactly at the point of the roughening transition, where 
large-scale fluctuations of diverge logarithmically and 
can be described by a (dimensionless) free-field effective 
Hamiltonian, 

= ^ y d'viVu^f , (A3) 

with dimensionless rigidity modulus J — tt/8. This value 
of J ensures the marginality of the operator — cos(7ru"^) 
favoring even values of w^. It follows then from Eq. 
(Ala) that the correlation function (c^crj^) decays al- 
gebraically, {o-f^crfj (X l/rl^ with r]^ = 1, ri2 being the 
distance between ji and j2- 

It is clear already from symmetry that the analogous 
correlation function on Tb has to decay in exactly the 
same way as on 7a . The law of this decay can also be 
derived with the help of another approach. It can be rig- 
orously shown that the correlation function of the two 
pseudospins on sublattice Tb is equal to the probability 
that they belong to the same BC loop. In terms of SOS 
representation such a loop corresponds to a closed step 
at which variables uf jumps by ±2. 

In the framework of the SOS model defined on sub- 
lattice Ta, the total statistical weight of configurations 
in which the two given sites on Tb belong to the same 
BC loop can be found by reversing the sign of the step 
on one of the two segments into which these two sites 
split the loop.^^'"*" This transforms the full set of such 
configurations into the set of configurations allowed in a 
system with a neutral pair of screw dislocations on go- 
ing around each of which changes hy b = ±4. When 
the large-scale fluctuations of u"^ can be described by 
the Hamiltonian (A3), the dimensionless free energy as- 
sociated with such a dislocation pair will be given by 
F « (J6^/27r) lnri2, which shows that (c^cr^) ocr^2- 



The same value of the exponent follows also from the pos- 
sibility to unambiguously interpret the steps in the con- 
sidered model as the contours of equal height.''^ The con- 
sistency of the results produced by different approaches 
supplies an additional confirmation that the value of the 
rigidity modulus J in Eq. (A3) has been correctly chosen. 

In the framework of the reduced problem, the ex- 
pression for the chirality of a plaquette is reduced to 
the product of two pseudospins on neighboring sites, 
Xk = '^jA{k)'^jB{k)- I'^ order to find how chirality corre- 
lations decay with distance, it is convenient to express 
this product in terms of locally coarse-grained height 
f^k = hi^jAik) +Ujs{k}] as 

Xk = Cj^{k)<Jj^(k) = sin(7r/ife) . (A4) 

The variables are defined on AB bonds and acquire 
values which are shifted from integers by 1/2. It follows 
from the definition of hk that the values of these variables 
on adjacent bonds can be either equal to each other or 
differ by ±1. In particular, in any regular state the val- 
ues of hk are the same on all AB bonds. Since the large- 
scale fluctuations of h have to be described by the same 
Hamiltonian, Eq. (A3), as those of u"^, it follows from Eq. 
(A4) that chirality correlations have to decay with expo- 
nent rj^ = A (as in the four-state antiferomagnetic Potts 
model, see Sec. IV). This means that x is a marginal 
operator. 

Note that the decay of {xkXk') is much faster than that 

of the product of {(y]A(k)(^3A(k')) and ((TjB(fc)^iB(fe'))- 
terms of the loop representation the origin of this prop- 
erty is quite clear. The presence of an AC loop passing 
through the sites jA{k) and jA{k') strongly decreases the 
number of configurations with a BC loop passing through 
the sites jB{k) and jB{k'). From the analogy with the 
case of a single loop (discussed two paragraphs above), 
the fraction of configurations in which two such loops are 
present simultaneously can be estimated by considering 
a system with a neutral pair of screw dislocations with 
b — ±8, which again gives 77^ = 4. 
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